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Abstract 

The best natural candidates for the realization of color superconductivity are quark stars -not 
yet confirmed by observation- and the extremely dense cores of compact stars, many of which have 
very large magnetic fields. To reliably predict astrophysical signatures of color superconductivity, 
a better understanding of the role of the star's magnetic field in the color superconducting phase 
that realizes in the core is required. This paper is an initial step in that direction. The field 
scales at which the different magnetic phases of a color superconductor with three quark flavors 
can be realized are investigated. Coming from weak to strong fields, the system undergoes first 
a symmetry transmutation from a Color-Flavor-Locked (CFL) phase to a Magnetic-CFL (MCFL) 
phase, and then a phase transition from the MCFL phase to the Paramagnetic-CFL (PCFL) phase. 
The low-energy effective theory for the excitations of the diquark condensate in the presence of a 
magnetic field is derived using a covariant representation that takes into account all the Lorentz 
structures contributing at low energy. The field-induced masses of the charged mesons and the 
threshold field at which the CFL — > MCFL symmetry transmutation occurs are obtained in the 
framework of this low-energy effective theory. The relevance of the different magnetic phases for 
the physics of compact stars is discussed. 

PACS numbers: 12.38.Aw, 24.85. +p. 26.60,+c 
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I. INTRODUCTION 



At present the physics community is actively trying to find ways to differentiate a neutron 
star made up entirely of nuclear matter from one with a color superconducting [l| quark 
matter core. Compact stars typically have very large magnetic fields. Hence any predicted 
signature of a color superconducting core should take into account the presence of the star's 
magnetic field and its effects in the superconducting state. Given that magnetars can have 
surface fields as large as 10 14 — 10 16 G |2j, it is reasonable to expect that the star's interior 
fields can reach even higher values. Maximum strengths of 10 18 — 10 19 G are allowed by a 
simple application of the virial theorem 

Although a color superconductor (CS) is in principle an electric superconductor, because 
;he diquark condensate carries nonzero electric charge, in the color-flavor-locked (CFL) phase 
4| (CFL is the color superconducting phase that realizes in a system of three-flavor massless 
quarks at high densities) there is no Meissner effect for a new in-medium electromagnetic 
field A^. This in-medium electromagnetic field -called in the literature a "rotated" electro- 
magnetic field, where the "rotation" takes place in an inner space- is a combination of the 
regular electromagnetic field and the 8 th gluon 4], |5| . As the quark pairs are all neutral with 
respect to the "rotated" electromagnetic charge Q, the "rotated" electromagnetic field 
remains long-range within the superconductor. 

In this paper we are interested in the color-superconducting magnetic phases that realize 
in a very dense system of three-flavor massless quarks interacting in the background of a 
rotated magnetic field B. As shown in Ref p], the color-superconducting properties of such 
a system are substantially affected by the penetrating B field and as a consequence, a new 
phase, called Magnetic Color Flavor Locked (MCFL) phase (J, takes place. In the MCFL 
phase the pairing of (rotated) electrically charged quarks is reinforced by the field. Pairs of 
this kind have bounding energies which depend on the magnetic-field strength and are bigger 
than the ones existing at zero field. At field strengths of the order of the baryon chemical 
potential, the pairing reinforcement is sufficient to produce a distinguishable splitting of the 
gap in two pieces: one that only gets contributions from pairs of neutral quarks and one 
that gets contributions from both pairs of neutral and pairs of charged quarks. 

Although the symmetry breaking patterns of the MCFL and CFL phases are different, 
the two phases are hardly distinguishable at weak magnetic fields. In the CFL phase the 
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symmetry breaking is given by 



G = SU{3) C x SU(3) L x SU(3) R x U(1) B x U(l) e . m . -> Sf/(3) c+L+R x C/(l) e . m .. (1) 

This symmetry reduction leaves nine Goldstone bosons: a singlet associated to the breaking 
of the baryonic symmetry U(1)b, and an octet associated to the axial SU(3)a group. 

Once a magnetic field is switched on, the difference between the electric charge of the u 
quark and that of the d and s quarks reduces the original flavor symmetry of the theory and 
consequently also the symmetry group remaining after the diquark condensate is formed. 



Then, the br ea ki n g pat t e r n for the MCPL-phase^ beccnes 

Gb = SU(3) c xSU(2) L xSU(2) R xU(l)^xU(l) B xU(l) e . m . -> SU(2) c + L+ RxU(l) e . m , (2) 

The group U(1)a (not to be confused with the usual anomaly U(1)a) is related to the 
current which is an anomaly-free linear combination of s, d, and u axial currents (7|. In 
this case only five Goldstone bosons remain. Three of them correspond to the breaking of 
SU(2)a, one to the breaking of ?7(1)^ , and one to the breaking of U(1) B . Thus, an applied 
magnetic field reduces the number of Goldstone bosons in the superconducting phase, from 
nine to five. 

The MCFL phase is not just characterized by a smaller number of Goldstone fields, but 
by the fact that all these bosons are neutral with respect to the rotated electric charge. 
Hence, no charged low-energy excitation can be produced in the MCFL phase. This effect 
can be relevant for the low energy physics of a color superconducting star's core and hence 
for its transport properties. In particular, the cooling of a compact star is determined by 
the particles with the lowest energy; so a star with a core of quark matter and sufficiently 
large magnetic field can have a distinctive cooling process. 

More recently, we have found that the magnetic field can also influence the gluon dynamics 
js]. At field strengths comparable to the charged gluon Meissner mass a new phase can be 
realized giving rise to an inhomogeneous condensate of Q-charged gluons [8|. The gluon 
condensate anti-screens the magnetic field due to the anomalous magnetic moment of these 
spin-1 particles. Because of the anti-screening, this condensate does not give a mass to the 
Q photon, but instead amplifies the applied rotated magnetic field. This means that at such 
applied fields the CS behaves as a paramagnet, thus we named this phase paramagnetic 
CFL (PCFL) This last effect is also of interest for astrophysics. Compact stars with 



color superconducting cores could have larger magnetic fields than neutron stars made up 
entirely of nuclear matter, thanks to the gluon vortex antiscreening mechanism. 

The above state of affairs underlines the need to discern the scales and field strengths 
at which one or another magnetic phase is physically relevant. If one ignores the quark 
masses, the main scales of the color superconductor are the baryon chemical potential fi, 
the dynamically generated gluon mass m g ~ g/i and the gap parameter A ~ ^e~ a ^ 9 , with 



a a constant that is dominated by magnetic gluon exchanges [10] . We can assume that at 
sufficiently high fi, the running strong coupling g becomes <C 1, so the the hierarchy of 
the scales is A <C m g <C /i. The main purpose of this paper is to elucidate how the different 
magnetic phases are related to the fundamental scales of the CS. 

The plan of the paper is as follows. In Sec. [Ill we develop the CFL low-energy effective 
theory in a magnetic background using a Lorentz covariant formalism. In this derivation the 
rotated magnetic field is introduced only through covariant derivatives, thus preserving the 
U{l) em gauge invariance. The threshold rotated magnetic field that decouples the charged 
mesons from the low-energy theory is found in Sec. IIHI In Sec. [IV] we discuss how the 
different magnetic phases are realized in a hierarchical order determined by the main energy 
scales of the CS. In the concluding remarks we state the major outcomes of the paper 
and discuss possible astrophysical implications of the realization of a PCFL-like phase at 
moderate densities. 



II. LOW-ENERGY EFFECTIVE THEORY IN A MAGNETIC BACKGROUND 

The physics at energies below the lower scale A can be explored by constructing the 
effective low-energy theory in the presence of a rotated magnetic field. Since in the CFL 
phase all the fermions are gapped and all the gluons have dynamically generated masses, the 
low energy theory of the CS is governed by the Goldstone modes arising from the breaking 
of the global symmetries in the presence of a magnetic field. 

As discussed in the Introduction, once a magnetic field is present, the original symmetry 
group Q is reduced, due to the different electric charges of the quarks, to Qg. One would 
think that the low energy theory should correspond to the breaking pattern ([2j), hence be 
described by five neutral Goldstone bosons. However, it is clear that at very weak magnetic 
fields the symmetry of the CFL phase can be treated as a good approximated symmetry, 
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meaning that at weak fields the low-energy excitations are essentially governed by nine 
approximately massless scalars (those of the breaking pattern ([T])) instead of five. 

A question of order here is: what do we exactly understand as a very weak magnetic 
field? In other words, what is the threshold-field strength that effectively separates the CFL 
low energy behavior from the MCFL one? A fundamental clue in this direction will come 
from the determination of the term in the low-energy CFL Lagrangian that can generate a 
field-induced mass for the charged Goldstone fields, disconnecting them from the low-energy 
dynamics at some field strength and thereby effectively reducing the number of Goldstone 
bosons from the nine of the CFL phase, to the five of the MCFL. A similar approach was 

n 

previously followed in Ref. Nevertheless, as it will become clear below, our treatment 
and results will differ from those previously obtained. 

Our strategy will consist of writing the effective low energy Lagrangian for the Goldstone 
bosons corresponding to the CFL breaking pattern (pQ), but in the presence of an external B 
field, that is, ignoring the explicit breaking introduced by the electromagnetic interaction. To 
ensure that this Lagrangian is "invariant" with respect to the original group of the CFL case, 
one treats the charge operator as a spurion field 11], and assume it transforms conveniently 
under left /right flavor symmetries, as well as under the color symmetries, depending on 
whether the operator appears with flavor or color indexes, as it will be shown below. 

There are two important points that separate our treatment from previous works. One 
is that we give the general form of the low-energy Lagrangian in an arbitrary reference 
frame, that is, we introduce all the possible covariant structures that can be formed at finite 
density and in the presence of an external magnetic field. When taken in the rest frame, 
the Lagrangian naturally reproduces the different Lorentz structures that characterize the 
problem with external magnetic field at finite density. The other is that when proposing 
the allowed terms of the low energy Lagrangian, we take into account that the coupling of 
the charged mesons with a rotated electromagnetic field can be traced back to the coupling 
of the fermions with the field. This coupling only occurs within a covariant derivative to 
preserve the gauge invariance of the theory under the U(l) group. Then, the coupling 
between the Goldstone bosons and the rotated electromagnetic field always occurs within 
a covariant derivative too. We will show that in this U(l) gauge invariant approach, the 
charged Goldstone bosons acquire field-induced masses that appear at the leading order of 
the low-energy theory. 
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To find the effective low-energy theory, we can follow a similar method to that used 

n n 

at zero field in Refs. [12|, H3J . We start by introducing two scalar fields which describe the 
fluctuations of the diquark condensates, and are associated to left and right order parameters: 

X ai _ e *bc e ijk^b^ck^ ^ yai _ ^ijk (^g^)t (3) 

where i,j, k denote flavor indices, a,b,c denote color indices, and L/R denote left/right 
chirality, respectively. Under an SU(3)c x SU(3)l x SU(3)r rotation, the above fields 
transform as 

X - g c Xgl , Y - g c Y g\ , (4) 

with gc G SU(3)c, and G SU(3)r/l- The expectation values of the X and Y fields 
define the ground state of the CFL phase which produces the symmetry breaking SU(3)c X 
SU(3) L x 5C/(3) fl -> Sf/(3) c+L+fl . 

We are interested in the fluctuations of the phases of the order parameters. Therefore, 
we can factor out the norm of the order parameters in (jlj) and work with unitary X and Y 
scalar fields. Although the axial group U(1)a is anomalous at low densities, it becomes an 
approximate symmetry at high densities due to the suppression of the instanton interactions 
that produce the anomaly. Nevertheless, in our derivations neither the pseudo Goldstone 
mode associated to the breaking of this group, nor the Goldstone mode associated to the 
[/(l)-baryon symmetry breaking will be considered, as these bosons are both neutral with 
respect to the rotated electromagnetic charge, and therefore irrelevant for the analysis of 
the present paper (we refer the interested reader to Ref. 13j to find out the contributions 



of these modes to the low energy theory; see also Ref. [14| for the diquark excitations of the 
CFL ground state at zero magnetic field in the framework of a Nambu-Jona-Lasinio model). 
We can introduce the Goldstone canonical fields 

X = e m x T ", Y = e m ^, a = l,...8 (5) 

with T a being the SU(3) generators normalized to satisfy 

Tr[T a ,T b ]= l -5 ab (6) 

Thus, X and Y together contain a total of sixteen scalar fields. Strictly speaking only 
eight of them are genuine Goldstone bosons, that is, massless scalar fields associated to the 
breaking of global symmetries. The other eight are Higgs fields related to the spontaneous 



breaking of the color gauge group SU(3) that gives mass to the gluons and therefore can be 



always eliminated from the theory by choosing a convenient (unitary) gauge [13]. Since these 
considerations remain valid in the presence of the external magnetic field, we will work in 
the unitary gauge and concentrate our analysis into the derivation of the low energy theory 
for the genuine Goldstone bosons. 

At zero magnetic field, the low energy theory of the Goldstone bosons associated to the 



global symmetries can be written 



12| as 



4 

with currents J x and J Y defined by 



c = - f -fTr[{.r x - Jyf] (7) 



J x = X^X 1 ", J Y = Yd»Y*. (8) 

Notice that (jTJ) does not take into account the breaking of the Lorentz invariance due to the 
finite density. The lack of Lorentz invariance was later incorporated in Ref. 12| by manually 
assigning different coefficients in front of the temporal and spatial derivative terms. 
Using (jlj), we can verify that the currents Jx and Jy transform as 

J x ~> 9cJx9c> J y ~> 9cJy9c ( 9 ) 

The Lagrangian density (JTj) contains the leading (second) order in derivatives terms in- 
variant under the SU(3)c X SU(3)l x SU(3)r rotations. 

Two important changes occur when a B field is switch on. First, the derivatives should 
be replaced by covariant derivatives containing the rotated electromagnetic potential 
associated to the magnetic field B. Second, the number of fundamental tensors available in 
the theory increases, because of the extra tensor F^ u . As a consequence, we can construct 
an effective theory containing a larger number of independent terms which are quadratic in 
the (covariant) derivative. 

The covariant derivative should be consistent with the fact that the VEV of the diquark 
fields, which only get contribution from the diagonal elements of the order parameter matrix, 
are all neutral in the rotated charge. We can generically denote the 3x3 order parameter 
matrix by A. Taking into account that a rotated electromagnetic field is a combination of 
the 8 t/l -gluon field and the conventional electromagnetic field and that each element of the 
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order parameter matrix carries an electric charge equal to the sum of the electric charges of 
the quarks forming the corresponding pair, we can write the covariant derivative as 

[dp + igGl(T 8 x 1) + ieA,(l x Q A )]A = [d, - l -^-gGl(Q x 1) - ieA^l x Q)]A (10) 

where the direct products denote (color x flavor). Q<\ is the conventional electromagnetic 
charge operator of the quark pairs in the quark representation (s, d, u). In this representation 
the first column of A will have (d, u) pairs, the second (s, u) pairs, and the third (s, d) pairs. 
Hence Qa = diag (1/3, 1/3, —2/3). Both the T 8 generator of the SU(3) group and the usual 
quark charge operator Q = diag(— 1/3, —1/3, 2/3) are connected to Qa through the relations 
T 8 = (V3/2)Q A = -( v / 3/2)Q. 

Taking into account that A^ = A^cosO — G^sinO , where 9 is the mixing angle 4|, Eq.( flO|) 
can be rewritten as 

(d^-ieA^A (11) 

where e = ecosO = (y/3/2)g sinO, and Q — (1 x Q — Q x 1) is the "rotated" electric 
charge operator. As expected, the charge operator Q assigns zero rotated charge to the 
diagonal elements of A. Notice that we do not keep the orthogonal combination G 8 = 
A^sinO + G^cos9, as the field G 8 acquires a mass larger than A C fl, so it decouples from 
the low-energy theory. 

The straightforward generalization of the currents (JHJ) to the case with nonzero rotated 
magnetic field should be done by substituting the derivative in (jHJ) by the covariant derivative 
( TTTj) . Taking into account that XX^ = 1 one can write the left current as 

( J£U = Xai^Sij + ieA^Q l3 )X] b - ieA»Q ab = X(d»X^ + ieA»Q L X^) - ieMQ c (12) 

where we introduced the notations Q c and Q L to keep track of whether the operator Q is 
an operator in color or (left) flavor space. A similar expression can be found for the right 
current Jy, with the obvious substitution X — > Y and Q L — > Q R . 

As mentioned above, once a magnetic field is present one has an extra tensor in the 
system that allows to create new structures in the Lorentz space. Working in a covariant 
way, the leading order low-energy Lagrangian density can be written as 

c = -S-Tr[(r x - 4)(r x - r Y )Q, v ] (is) 



s 



with 

©/^ = Ci5V + C 2 u fjL u u + C 3 F pp F pu (14) 

being the most general Lorentz structure that can be formed with the vectors and tensors 
available at low energies. In (1141) we considered the normalized electromagnetic tensor 
Ffj, v = rmF^, in addition to the usual metric tensor g pv and the vector four-velocity of 
the center of mass of the dense medium u p . The C[s are just constant coefficients. In 
this covariant representation the magnetic field can be expressed as B p = \e pi/p \u v F pX . In 
the rest frame, u u = (1,0,0,0), and the electromagnetic tensor becomes F pu = BS^id^ if 
we assume a magnetic field pointing along the third spatial direction. In the absence of 
a magnetic field, the coefficient C 3 is taken equal to zero and the Lagrangian density ( |T3l) 
reduces in the rest frame to that introduced in Refs. {12I . Il3j |. where the system was not 
Lorentz invariant due to the finite density (u p 7^ 0), but it kept the rotational symmetry. In 
the presence of a magnetic field the structure associated to C3 naturally separates between 
modes longitudinal and transverse to the field. A linear term in F pu is forbidden, since 
it would violate the theory CP invariance. Notice that in the leading (second) order in 
derivatives we do not need to introduce the momentum k p in the structures contributing to 

(ED. 

Notice that the flavor symmetry SU(3)l x SU(3)r is explicitly broken by the electro- 
magnetic coupling in (|T3|) . So, strictly speaking, (jT3l) is invariant under Q B but not under 
Q. However, one can make the theory invariant under the original group Q if we treat the 
charge operators as spurion fields and assume that they transform as 

Q L - 9lQ L gl Q R - g* R Q R gl, Q c - 9cQ c gl (15) 

Using these transformations, one can show that the new currents transform under Q in the 
same way as the currents at zero field: 

Jx ~* 9cJx9c, Jy -> 9cJy9c ( 16 ) 

yielding to the invariance of (|T3j) under Q. 

Following Ref. [3], we introduce now the color singlet, 

£ = r f x, (17) 
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which transforms under SU(3)c x SU(3)l x SU(3)r as E ^ g* F ^lgJ j . In terms of S the 
Lagrangian density f|T3|) can be written as 

C = f -^Tr[{D^){D^Q pv \ (18) 

where the covariant derivative acting on X is 

= <9 M £ + ieA M (Q fl £ - SQ L ) (19) 

In the rest-frame, ffT8l) becomes 



£ = ^ [Tr (D Z)(D tf) +v 2 ± Ti (D ± X)(D ± tf) + ujTr (£>„£)p[|£t)] . (20) 

showing the expected separation between longitudinal Dy and transverse Di (to the field) 
components of the covariant derivatives, a direct consequence of the partial breaking of the 
rotational symmetry in the presence of the magnetic field. The decay constant /„., and the 
meson maximum velocities v±, and v\\, are parameters to be computed from the microscopic 
theory. They can in general depend on the baryonic chemical potential and the applied 
magnetic field. In the weak-field limit (eB <C /i 2 ), at asymptotically large values of the 



chemical potential /a, they can be approximated by their zero- field values 13] 



111 ] (Eq. (2.16)), is present 



, 2 21-81n2 u 2 1 

s ^^^h> ^""i"^ < 21) 

As we shall see in the next Section, the low-energy theory (1201) supports the generation of 
field- dependent masses for the charged mesons at fields larger than a threshold value, even 
though no quadratic-in-5 term, like the one proposed in Ref. 
at the leading order. To understand this, one should keep in mind that the interaction 
between £ and the electromagnetic field always occur through the covariant derivative. To 
generate a quadratic-in-5 (and quadratic-in-Q) term, one needs to consider a higher order 
contribution involving four covariant derivatives 

J d 4 xTri(D^)(D u Z)\D p Z)(D x ZyF^F pX \ (22) 

Using (we dropped total derivatives) 

D^DjSfF^ = X -F llv \mF llv {Q R -ZQ L tf) + 2ieA^Q R Zd v tf -Zd„tfQ R )-2E(d u E^)d^ 

(23) 
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it can be straightforwardly shown that the term proposed in 

J d 4 x?B 2 Tr[Q R ^Q L ^} (24) 

is just one of the several terms coming out of (|22|) . However, as proved below, the threshold 
field for the decoupling of the charged mesons will be determined by a contribution more 
relevant than AMD- 



III. CFL-MCFL THRESHOLD FIELD 



The unitary matrix X can be parametrized in term of the elementary Goldstone bosons 



y 4 as 



l^ArpA 



E = expU^^J , A=l,...,8, (25) 

where T A are the £77(3) generators. Expanding (1251) up to linear terms in the fields, we can 
write (fT9|) as a sum of covariant derivatives for the (f) A fields 

3 

= f E tAq ^ a + T%d ^ + E r± ^ ± ^AJ n± + E A± ^ ± ^) K± l ( 26 ) 

^ A=l ± ± 

where 

Il ± = ^Tici>% K+^^^if), (27) 

r ± = 4= ± , A ± = [T 6 ± *T 7 ] (28) 
v2 v2 

and the rotated electromagnetic potential is taken in the Landau gauge 

Au — (0, 0, Bx, 0) (29) 
From (TSUI) , the low-energy Lagrangian for the Goldstone bosons can be written as 

L = I ^E W A I 2 + l 9 ° T H + f E I^i^T + l 5 n T l 2 ] 

J A=\ A=l 

2 3,8 

+TE \ d ^\ 2 + + ^±)T| 2 ]} (30) 

A=l 

where we introduced the charged meson doublet 
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Y = , (31) 



The Lagrangian (1301) represents the CFL low-energy theory in the presence of a weak 
(k 2 ~ eB <C /i 2 ) constant magnetic field. Now we are ready to determine the strength of the 
threshold field for which the effective symmetry transmutation from CFL to MCFL occurs. 

For that aim, it is convenient to work in momentum space. Transforming to momentum 
space in the presence of the magnetic field can be done by applying to the scalar fields the 
same method originally developed for fermions in [15| and later extended to vector fields in 
In this approach the transformation to momentum space can be carried out using the 
wave functions Sk{x) of the asymptotic states of the charged mesons in a uniform magnetic 
field. These functions play the role in the magnetized medium of the usual plane-wave 
(Fourier) functions e tpx at zero field. Then, for the charged field T we have 



T(Jfe) = J d A xS k (x)T(x) (32) 

where 

Sk{x) — J\fexp(ik x° + ik 2 x 2 + ik 3 x 3 )D n (p) (33) 



with D n (p) being the parabolic cylinder functions with argument g = y 2\eB\(xi — k 2 /eB), 
J\f the normalization constant (J\f = (47r|eB|)*/\/nF), and n = 0, 1, ... denoting the Landau 
levels. It is easy to check that the transformation functions (l33jl satisfy the orthonormality 
condition 

/CO 
d 4 xS k (x)S k ,(x) = (27r) 4 ?(fc / - k) (34) 
-00 

where 8{k' — k) = 5(k' — k )5(k 2 — k 2 )5(k' 3 — k 3 )5 n ' n . 

Using this transformation in ( 1301) we can derive the Klein-Gordon equation for the charged 
mesons in momentum space 

[k 2 - eB{2n + l)v\ - klv 2 }T(k) = 0. (35) 
and from it the corresponding dispersion relation 

E 2 = eB{2n + l)v\ + k 2 3 vl (36) 
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We see that at zero momentum (k^ = 0, n = 0) the rest energy of the charged mesons is 
M~ = eBvj_, meaning they acquire a mass induced by the magnetic field. For a meson to be 
stable, its mass should be less than twice the gap, otherwise it will decay into a quasiparticle- 
quasihole pair. Here we are assuming that the interactions between the Goldstone bosons and 
the quasiparticles are described by a Yukawa term that originates within a NJL-type model. 
Son, authors , S ee Q for deta*) have ar gU ed that the ndcroscopic structure of the NG 
bosons in the CFL phase should be that of a quartic-quark state. In the quartic-quark picture 
the threshold would presumably be four times the gap energy. Within the NJL-approach, 
the threshold field where the number of stable Goldstone bosons effectively changes, thus the 
symmetry of the CFL phase can not longer been treated as a good approximated symmetry 
is given by 

eB MCF L = \^cfl - 12A^ FL , (37) 

v ± 

where we used the weak-field approximation v± ~ l/\/3 jTsj] . Contrary to the result found in 
our threshold field does not depend on the decay constant f^, therefore it depends on fi 
only through Acfl- The f n dependence found in [HI] is a direct consequence of considering 
a subleading contribution in the derivation of the threshold field value, as previously shown. 

For Aqfl ~ 15MeV we get eBucFL ~ 10 16 G. At these field strengths, the charged 
mesons decouple from the low-energy theory. When this decoupling occurs, the five neutral 
Goldstone bosons (including the one associated to the baryon symmetry breaking) that 
characterize the MCFL phase will drive the low-energy physics of the system. Therefore, 
coming from low to higher fields, the first magnetic phase that effectively shows up in the 
magnetized system will be the MCFL, even though at fields near the threshold field the 
splitting of the gaps found at much stronger fields [6| may still be negligible. 

We underline that the phenomenon occurring at the threshold field (1371) is not a phase 
transition, as no symmetry is broken there. At any nonzero magnetic field strength, below 
or above the threshold field ( |371) . the symmetry of the system is, strictly speaking, that of 
the MCFL. However, at B < Bmcfl the charged Goldstone bosons are so light that the 
observable impact of the smaller symmetry of the MCFL phase, compared to that of CFL, 
is irrelevant and the nine Goldstone bosons of the CFL phase provide a good approximated 
description of the low-energy physics. Based on these considerations we choose to call the 
CFL-MCFL transition a symmetry transmutation. 
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On the other hand, it is worth to call attention to the analogy between the CFL-MCFL 
transmutation and what could be called a "field- induced" Mott transition. Mott transitions 
were originally considered in condensed matter in the context of metal-insulator transitions 
in strongly-correlated systems [la ]. Later on, Mott transitions have been also discussed in 
QCD to describe derealization of bound states into their constituents at a temperature 



defined as the Mott temperature [19|. By definition, the Mott temperature Tm is the tem- 
perature at which the mass of the bound state equals the mass of its constituents, so the 
bound state becomes a resonance at T > Tm- In the present work, the role of the Mott 
temperature is played by the threshold field Bmcfl- Mott transitions typically lead to the 
appearance of singularities at T = Tm in a number of physically relevant observables. It is 
an open question worth to be investigated whether similar singularities are or not present 
in the CFL-MCFL transmutation at Bmcfl- 

IV. B VS fi PHASES IN A COLOR SUPERCONDUCTOR WITH THREE QUARK 
FLAVORS 

What will happen if we keep increasing the magnetic field until it reaches the next energy 
scale g/il As known [8], due to the interaction of the applied magnetic field with the charged 
gluon anomalous magnetic moment (ief^ u G+G~), once B > B PC fl = m M> with m M ~ gf-i 
being the magnetic mass of the charged gluons, one of the modes of the charged gauge field 
becomes tachyonic (this is the well known "zero- mode problem" for spin-1 charged fields in 



the presence of a magnetic field found for Yang-Mills fields 20], for the bosons in the 



electroweak theory 



21 



22], and even for higher-spin fields in the context of string theory 



23l|) . Similarly to other spin-1 theories with magnetic instabilities 20]]- 22J, the solution of 
the zero-mode problem leads to the restructuring of the ground state through the formation 
of an inhomogeneous gauge-field condensate G, as well as an induced magnetic field due to 
the back reaction of the G condensate on the rotated electromagnetic field. The magnitude 
of the G-condensate plays the role of the order parameter for the phase transition occurring 
at B = Bpcfl- Near the transition point, the amplitude of the condensate G is very small 



,8 



. Then, the condensate solution can be found using a Ginzburg-Landau (GL) approach 
similar to Abrikosov's treatment of type II metal superconductivity near the critical field 
H c2 [24]. As in Abrikosov's case, the order parameter \G\ continuously increases from zero 
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with the applied magnetic field, signalizing a second-order phase transition towards a gluon 
crystalline vortex state characterized by the formation of flux tubes. Both spatial symmetries 
-the rotational symmetry in the plane perpendicular to the applied magnetic field and the 
translational symmetry- are broken by the vortex state. 

It turns out that, contrary to what occurs in conventional type-II superconductors, where 
the applied magnetic field only penetrates through flux tubes and with a smaller strength 
than that of the applied field, the gluon vortex state exhibits a paramagnetic behavior. That 
is, outside the flux tube the applied field B totally penetrates the sample, while inside the 
tubes the magnetic field becomes larger than B. This antiscreening behavior is similar to 
that found in the electroweak system at high magnetic field 22j. Hence, since the Q photons 
remain long-range in the presence of the condensate G, the U(l) em symmetry remains un- 
broken. At asymptotically large densities, because Acfl <C tum, we have Bmcfl <C Bpcfl 
for each \i value. 

At fields B > fi 2 the density of states on the Fermi surface of the charged quarks will 
be larger than that of neutral quarks 6]. Because of the different density of states, the 
magnitude of the gap receiving contributions from pairs of rotated charged quarks will split 
from the magnitude of the gap receiving contributions only from pairs of rotated neutral 
quarks, as shown to happen, without taking into account the vortex state, in Ref. jf]]. The 
splitting of the gaps at this scale, however, would not break any new symmetry that has not 
been already broken at much lower scales. Hence, no new phase transition occurs at fields 
of order fi 2 . 

For fields much larger than fi 2 , i.e. sufficiently strong as to surpass all the energy scales of 
the system, the quark infrared dynamics will become predominant, and the phenomenon of 



magnetic catalysis of chiral symmetry breaking |25| will be activated, producing a phase that 
favors quark-antiquark condensates over quark-quark condensates. However, the exploration 
of this region goes beyond the scope of this paper. 

FigJT] provides a qualitative sketch of the magnetic phases that exist at asymptotically 
high densities in the framework of a three-flavor color superconductor. In that region, at 
very weak magnetic fields, the color superconducting state is practically described by the 
CFL phase, because the charged mesons, although massive, are so light that they cannot 
decay in pairs of quasiparticle-quasihole. When the field strength is of the order of the 
quarks' energy gap, the charged mesons become heavy enough to decouple and the low- 
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FIG. 1: Qualitative sketch in the B vs fi 2 plane of the different phases of a color superconductor with 
three quark flavors in the presence of an external magnetic field at asymptotically high densities. 
The CFL phase appears here as an approximate symmetry at weak field. Thus, the line between 
the CFL and MCFL phases does not denote a real phase transition, but the boundary separating 
the approximated CFL phase from the MCFL phase. This symmetry-transmutation line is reached 
at field values of the order of the CFL gap square. The line between the MCFL and PCFL phases 
indicates a second-order phase transition curve occurring at field strengths of the order of the 
magnetic mass square of the charged gluons. The rectangular region to the left corresponds to 
moderately high densities in the presence of a magnetic field. Since the ground state at moderately 
high densities has not yet been investigated in the presence of a magnetic field, this region is 
indicated by question marks. 



energy physics is indeed that of the MCFL phase, where five neutral massless bosons drive 
the low-energy behavior. At fields comparable to the magnetic masses of the charged gluons, 
a chromomagnetic instability is developed for these gluons leading to the formation of a 
vortex state and the antiscreening of the magnetic field The vortex state breaks 

the translational symmetry, as well as the remaining rotational symmetry in the plane 
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perpendicular to the applied magnetic field, hence the vortex formation corresponds to a 
phase transition from the MCFL to a PCFL phase. 

The chemical potential Hcfl in Fig{T]is used to schematically separate the regions where 
the effects of an s-quark mass M s can (right of hcfl) ot cannot (left of Hcfl) be neglected. 
The region of moderately high densities to the left of [Iqfl alon g th e zero-magnetic- field 
line has been subjected to intense scrutiny in the literature 26j]-32|], as this is the most 
important region for applications of CS in realistic systems as compact stars. A main 
problem has been to find the stable phase at these moderately high densities. Despite many 
clever propositions that include a modified CFL-phase with a condensate of kaons 29j; 
a LOFF phase on which the quarks pair with nonzero total momentum [30]; as well as 



homogeneous 



31] and inhomogeneous [32] gluon condensate phases, it is still unclear which 



of these phases produces the lowest free energy. Given that nobody has yet studied this 
question at finite magnetic field, we choose to indicate it in the figure with question marks. 



V. CONCLUDING REMARKS 



Summarizing, in a color superconductor with three-flavor quarks at very high densities 
an increasing magnetic field produces a phase transmutation from CFL to MCFL first, 
and then a phase transition from MCFL to PCFL. During the phase transmutation no 
symmetry breaking occurs, since in principle once a magnetic field is present the symmetry is 
theoretically that of the MCFL, as discussed above. However, in practice for B < B MC fl ~ 
A"cfl the MCFL phase is almost indistinguishable from the CFL. Only at fields comparable 
to ^cfl the main features of MCFL emerge through the low-energy behavior of the system. 
At the threshold field Bmcfl, only five neutral Goldstone bosons remain out of the original 
nine characterizing the low-energy behavior of the CFL phase. These are precisely the 
five Goldstone bosons determining the new low-energy behavior of the genuinely realized 
MCFL phase. Going from MCFL to PCFL is, on the other hand, a real phase transition Q], 
as the translational symmetry, as well as the remaining rotational symmetry in the plane 
perpendicular to the applied magnetic field are broken by the vortex state. 

Throughout this paper we have ignored the effects due to quark masses because we as- 
sumed very large baryon density. However, the densities of interest for most astrophysical 
applications are just moderately high. At moderate densities the s-quark mass can and do 
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play an important role 26] . In this case, the color superconductor develops chromomagnetic 
instabilities even in the absence of an external magnetic field 27|. At even lower densities, 
where the s-quark decouples due to its larger mass, a two-flavor color superconductivity -the 
so-called 2CS phase- realizes. In this phase, when color neutrality and (3 equilibrium con- 
ditions are imposed, some chromomagnetic instabilities can also develop at certain density 
values [28(. Finding the stable superconducting ground state at moderate densities is one 



of the main questions in the field at the present moment [29(- 32|. The understanding of 
this problem in the presence of a magnetic field, which no doubt is another crucial player 
in the physical scenario of a compact star, is also an important open question. In this re- 
gard, as discussed in Ref. 32j, the removal of the chromomagnetic instabilities found at 
moderate densities in the two-flavor system may be related to the spontaneous generation 
of an inhomogeneous gluon condensate with a corresponding induced magnetic field. If this 
proposition is proved to be correct also for the three-flavor case, the star's core could be in 
a PCFL-like magnetic phase, even if the core's original magnetic field is zero or relatively 
low. That is, the PCFL-like phase will not be triggered by instabilities produced at a critical 
value of some pre-existing inner magnetic field, but by instabilities connected to the interplay 
of the neutrality conditions and the s-quark mass at some baryon density. An interesting 
consequence of a PCFL core is that a star's core in this phase can generate and/or boost its 
inner magnetic field. Exploring the magnetic phases that realize at realistic densities is an 
important pending task. 

It is plausible that if compact stars are the natural playground for color superconductivity, 
the magnetic phases described in this paper, or more precisely, the version of these phases 
at more realistic densities, may be relevant for the physics of the core of highly magnetized 
compact objects like magnetars 0, [3], and the so-called Central Compact Objects (CCO) 



34j . CCO are point-like sources located near the center of supernova remnants that cannot 



be identified as active radio pulsars or magnetars [3J|. Some of them may have magnetar's 
strength B fields and much smaller radii. For typical neutron star masses ~ 1.4M , a smaller 
radii means a denser star's core; the denser the core, the greater the chance it can be in a 
color superconducting phase. Magnetars' surface magnetic fields are typically in the range 
of 10 14 G — 10 16 G, and the fields at their much denser cores can probably reach even larger 
values. Any of these compact objects could be a candidate for the realization of a color 
superconducting phase in a highly magnetized background. 
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Moreover, the standard explanation of the origin of the magnetars' large magnetic fields 



cannot explain all the features of the supernova remnants surrounding these objects [351. 136] . 
Magnetars are supposed to be created by a magnetohydrodynamic-dynamo mechanism that 
amplifies a seed magnetic field due to a rapid rotating (spin period < 3ms) protoneutron 
star. Part of this rotating energy is supposed to power the supernova through rapid magnetic 
braking, implying that the supernova remnants associated with magnetars should be an order 
of magnitude more energetic than typical supernova remnants. However recent calculations 



35| indicate that their energies are similar. In addition, one would expect that when a 



magnetar spins down, the rotational energy output should go into a magnetized particle wind 
of ultrarelativistic electrons and positrons that radiate via synchrotron emission across the 
electromagnetic spectrum. Nevertheless, so far nobody has detected the expected luminous 
pulsar wind nebulae around magnetars 37]. 



Although more observations are needed to confirm the above, current observations indi- 
cate that alternative models to the standard magnetar model [2J need to be considered. For 
example, some authors [35] have suggested that magnetars could be the outcome of a stellar 
progenitor with highly magnetized cores. A progenitor star with a PCFL-like core would be 
capable to induce and/or enhance the star's magnetic field due to the antiscreening mecha- 
nism inherent to this color superconducting phase , and as such provides an alternative 
to the observational conundrum of the standard magnetar's paradigm 

Only after the stable color-superconducting phase at moderate densities is well established 
with and without an external magnetic field, we will be in a well-grounded position to 
investigate and reliably predict observable signatures of color superconductivity in compact 
stars. 
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